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Abstract 

In the framework of the Lindblad theory for open quantum systems, we determine 
the degree of quantum decoherence of a harmonic oscillator interacting with a thermal 
bath. It is found that the system manifests a quantum decoherence which is more and 
more significant in time. We calculate also the decoherence time scale and analyze the 
transition from quantum to classical behaviour of the considered system. 

1 Introduction 

The quantum to classical transition and classicality of quantum systems continue to be 
among the most interesting problems in many fields of physics, for both conceptual and 
experimental reasons [1, 2]. Two conditions are essential for the classicality of a quantum 
system [3]: a) quantum decoherence (QD), that means the irreversible, uncontrollable and 
persistent formation of a quantum correlation (entanglement) of the system with its envi- 
ronment [4], expressed by the damping of the coherences present in the quantum state of the 
system, when the off-diagonal elements of the density matrix decay below a certain level, 
so that this density matrix becomes approximately diagonal and b) classical correlations, 
expressed by the fact that the quantum state becomes peaked along a classical trajectory. 
Classicality is an emergent property of open quantum systems, since both main features of 
this process - QD and classical correlations - strongly depend on the interaction between the 
system and its external environment [1, 2]. In this work we study QD and analyze quantum- 
classical transition of a harmonic oscillator interacting with an environment, in particular 
with a thermal bath, in the framework of the Lindblad theory for open quantum systems. 



2 Master equation and density matrix 

In the Lindblad axiomatic formalism based on quantum dynamical semigroups, the irre- 
versible time evolution of an open system is described by the following general quantum 
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Markovian master equation for the density operator p(t) [5]: 

3 

The harmonic oscillator Hamiltonian H is chosen of the general quadratic form 

H = H + ^ qP + Pq ), H = ^-p* + ^f q z (2) 

and the operators Vj, Vj , which model the environment, are taken as linear polynomials in 
coordinate q and momentum p. Then the master equation (1) takes the following form [6]: 

dp % % % 

-j t = —fi[H ,p] - ^(A + fi)[q, pp + pp] + — (A - n)\p, pq + qp] 

-ff [?, k,p}} - ^[P, hp]] + ^([?, + [P, (3) 

The diffusion coefficients -D pp , D pg and the dissipation constant A satisfy the fundamental 
constraints: D pp > 0, D qq > and D pp D qq — D 2 q > \ 2 h 2 /4. In the particular case when the 
asymptotic state is a Gibbs state pc(oo) = e~~£r /Tre"fcr ; these coefficients become 



A + u. Tiu) \ — ah hu 

D m — Timuo coth — — -, D aa — coth — — , D va = 0, (4) 

pp 2 2kT' qq 2 mu 2kT' 



where T is the temperature of the thermal bath. In this case, the fundamental constraints 
are satisfied only if A > p and 

(A 2 - P 2 ) coth 2 ^ > A 2 (5) 

and the asymptotic values a qq (oo), a pp (oo), a pq (oo) of the dispersion (variance), respectively 
correlation (covariance), of the coordinate and momentum, reduce to [6] 

Oaa(oo) = coth——, <7 m} (oo) = coth——, er 00 (oo) = 0. (6) 

qqK ' 2mu 2kT' ppK ' 2 2kT' pqK ' K ' 

We consider a harmonic oscillator with an initial Gaussian wave function (<r q (0) and cr p (0) 
are the initial averaged position and momentum of the wave packet) 

m = ( ^(o) )} exp| -3^(5) (1 - !^ (0))( « - CT « (0))2 + l°> m < (7) 

representing a correlated coherent state (squeezed coherent states) with the variances and 
covariance of coordinate and momentum 

H5 hmu hr 

^ (0) = w appi0) = 2w^y apM = 2-vf^- (8) 
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5 is the squeezing parameter which measures the spread in the initial Gaussian packet and 
r, with |r| < 1 is the correlation coefficient. The initial values (8) correspond to a minimum 
uncertainty state, since they fulfil the generalized uncertainty relation a qq (0)a pp (0) — <7 pq (0) = 
h 2 /A. For 5=1 and r = the correlated coherent state becomes a Glauber coherent state. 
From Eq. (3) we derive the evolution equation in coordinate representation: 

dp ih d 2 d 2 imuj 2 2 /2 
dt = ^W~dq i2 ^^ ~ 

-i(A + ,.)(, - ,')<! - |> + i(A - „)[(, + 9 ')(| + A) + 2]p 
-^(9 - ?') 2 /> + O m (| + - 2iD r „fi(, - ,')(| + (9) 

The first two terms on the right-hand side of this equation generate the usual Liouvillian 
unitary evolution. The third and forth terms are the dissipative terms and have a damping 
effect (exchange of energy with environment). The last three are noise (diffusive) terms 
and produce fluctuation effects in the evolution of the system. D pp promotes diffusion in 
momentum and generates decoherence in coordinate q - it reduces the off-diagonal terms, 
responsible for correlations between spatially separated pieces of the wave packet. Similarly 
D qq promotes diffusion in coordinate and generates decoherence in momentum p. The D pq 
term is the so-called "anomalous diffusion" term and it does not generate decoherence. 
The density matrix solution of Eq. (9) has the general Gaussian form 

-W^) « - ^ + ~ ^ ~ 9<) + 5"M<« " ^ (10) 

where a(t) = a qq (t)a pp (t) — (J pq {t) is the Schrodinger generalized uncertainty function. In the 
case of a thermal bath we obtain the following steady state solution for t — > oo (e = ^): 

<^>-( S ^)W-^ + (.-^ a -.* (U, 
3 Decoherence and quantum- classical transition 

An isolated system has an unitary evolution and the coherence of the state is not lost - pure 
states evolve in time only to pure states. The QD phenomenon, that is the loss of coherence 
or the destruction of off-diagonal elements representing coherences between quantum states 
in the density matrix, can be achieved by introducing an interaction between the system and 
environment: an initial pure state with a density matrix which contains nonzero off-diagonal 
terms can non-unitarily evolve into a final mixed state with a diagonal density matrix. 
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Using new variables £ = (q + q')/2 and A = q — q', the density matrix (10) becomes 



p(E, A, i) = J- exp[-a£ 2 - 7A 2 + i(3Y,A + 2aa g (t)£ + - P<r q (t))A - aa 2 (t)}, (12) 

V 7i n 



with the abbreviations 



<r(t) Q _ a pq (t) 



2cr OT (t)' 7 ~ 2h 2 a qq (tY ha qq {ty 

The representation-independent measure of the degree of QD [3] is given by the ratio of 
the dispersion \/y/2rj of the off-diagonal element p(0, A,t) to the dispersion ^2/ 'a of the 
diagonal element p(£,0,t) : 

W*) = ^/?=-4fr- (14) 
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The finite temperature Schrodinger generalized uncertainty function has the expression 
[7] (with the notation Q 2 = u 2 — fi 2 , u > fi) 

h 2 1 

a(t) = — {e~ 4Xt [l- (S + — -) coth e + coth 2 e] 

4 5(1 -r 2 ) 

9X . , ... 1 n , uj 2 — a 2 cos(2Qt) 

+e~ 2Xt coth e [(5 + — ^ - 2 coth e A L 

0(1 — r z ) \V 

1 /xsnngQQ 2r/i^(l-cos(2^)) 1 2 

+(5 " W^) + nvi-^ ] + coth e} - (15) 

In the limit of long times Eq. (15) yields a (00) = (h 2 coth 2 e)/4, so that we obtain 

5qd(oo) = tanh^, (16) 

which for high T becomes <5qd(oo) = hu/(2kT). We see that 5qd decreases, and there- 
fore QD increases, with time and temperature, i.e. the density matrix becomes more and 
more diagonal at higher T and the contributions of the off-diagonal elements get smaller 
and smaller. At the same time the degree of purity decreases and the degree of mixedness 
increases with T. For T = the asymptotic (final) state is pure and 5qd reaches its initial 
maximum value 1. Sqd = when the quantum coherence is completely lost, and if 5qd — 1 
there is no QD. For long enough time the magnitude of the elements of the density matrix in 
the position basis are peaked preferentially along the diagonal q = q' . Then 5qd < 1 and we 
can say that the considered system interacting with the thermal bath manifests QD. Dissipa- 
tion promotes quantum coherences, whereas fluctuation (diffusion) reduces coherences and 
promotes QD. The balance of dissipation and fluctuation determines the final equilibrium 
value of 5qd- The initial pure state evolves approximately following the classical trajectory 
in phase space and becomes a quantum mixed state during the irreversible process of QD. 
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In the macroscopic limit, when H is small compared to other quantities with dimensions of 
action, the term in Eq. (9) containing D pp /h 2 dominates and induces the following evolution: 

| = -£fo-rtV <") 

Thus the density matrix loses off-diagonal terms in position representation, while the di- 
agonal (q = q') ones remain untouched. Quantum coherences decay exponentially and the 
decoherence time scale is of the order of 

h 2 

UeCO ~ D pp {q-q') 2 ' (18) 
In the case of a thermal bath, we obtain (see Eq. (4)) 

tdeco = 2 1 , (19) 

(A + ii)mwa gq (0) cothe' 

where we have taken (q — q') 2 of the order of the initial dispersion in coordinate cr qq (0). 
As expected, the decoherence time tdeco has the same scale as the time after which thermal 
fluctuations become comparable with quantum fluctuations. In the macroscopic domain QD 
occurs very much faster than relaxation. When t ^> t re i, where t re i ~ A -1 is the relax- 
ation time, which governs the rate of dissipation, the particle reaches equilibrium with the 
environment. Indeed, the uncertainty function a(t) (15) approaches a BE = (h 2 coth 2 e)/4, 
which is the Bose-Einstein relation for a system of bosons in equilibrium at temperature T. 
In the case of T = we approach the limit of pure quantum fluctuations, o"o = h 2 /A, which 
is the quantum Heisenberg relation. At high temperatures T we obtain the limit of pure 
thermal fluctuations, a MB = (kT/uo) 2 , which is a Maxwell-Boltzmann distribution for a sys- 
tem approaching a classical limit. The regime where thermal fluctuations begin to surpass 
quantum fluctuations is regarded as the transition point from quantum to classical statistical 
mechanics and the high temperature regime of a system is considered as the classical regime. 
We have shown that these two criteria of classicality are equivalent: the time when the 
quantum system decoheres is comparable with the time when thermal fluctuations overtake 
quantum fluctuations. After the decoherence time, the system has to be described by non- 
equilibrium quantum statistical mechanics. After the relaxation time the system is treated 
by equilibrium quantum statistical mechanics, and only at a sufficiently high temperature, 
when the spin statistics is represented by the Maxwell-Boltzmann distribution function, it 
can be considered in a classical regime [7, 8] . 

4 Summary and concluding remarks 

We have studied QD in the framework of the theory of open quantum systems in order 
to understand the quantum-classical transition for a harmonic oscillator in interaction with 
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a thermal bath. The classicality is conditioned by QD, expressed by the loss of quantum 
coherence in the case of a thermal bath at finite temperature. 

The role of QD became relevant in many interesting physical problems. In many cases one 
is interested in understanding QD because one wants to prevent decoherence from damaging 
quantum states and to protect the information stored in quantum states from the degrading 
effect of the interaction with the environment. Decoherence is also responsible for washing 
out the quantum interference effects which are desirable to be seen as signals in experiments. 
QD has a negative influence on many areas relying upon quantum coherence effects, in 
particular it is a major problem in the physics of quantum information and computation. 
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